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AGN Unification at 2 ~ 1: u — R Colors and Gradients in X-ray AGN Hosts 
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Magee^ 

ABSTRACT 

We present uncontaminated rest-frame u ~ R colors of 78 X-ray-selected AGN hosts at 
0.5 < z < 1.5 in the Chandra Deep Fields measured with HST/ACS/NICMOS and VLT/ISAAC 
imaging. We also present spatially-resolved NUV — R color gradients for a subsample of AGN 
hosts imaged by HST/WFC3. Integrated, uncorrected photometry is not reliable for comparing 
the mean properties of soft and hard AGN host galaxies at z ~ 1 due to color contamination 
from point-source AGN emission. We use a cloning simulation to develop a calibration between 
concentration and this color contamination and use this to correct host galaxy colors. 

The mean u—R color of the unobscured/soft hosts beyond ^ 6 kpc is statistically equivalent to 



o 
u 

rS^ . that of the obscured/hard hosts (the soft sources are 0.09 ±0.16 magnitudes bluer). Furthermore, 



the rest-frame V — J colors of the obscured and unobscured hosts beyond ~ 6 kpc are statistically 



O ' equivalent, suggesting that the two populations have similar distributions of dust extinction. For 

-*— > , the WFC3/IR sample, the mean NUV ~ R color gradients of unobscured and obscured sources 

j^ I differ by less than ~ 0.5 magnitudes for r > 1.1 kpc. These three observations imply that AGN 

obscuration is uncorrelated with the star formation rate beyond ~ 1 kpc. 

These observations favor a unification scenario for intermediate-luminosity AGNs in which ob- 
^ ' scuration is determined geometrically. Scenarios in which the majority of intermediate-luminosity 

f_^ ■ AGN at 2; ~ 1 are undergoing rapid, galaxy-wide quenching due to AGN-driven feedback pro- 

,__! i cesses are disfavored. 

(N 
l> 

O ; 1. INTRODUCTION 

The physical connection between embedded supermassive black holes and their host galaxies remains 
. ^^ _ an active topic of research. The disc overy of correlations between bulge properties and black hole mass (e.g. , 

j^ ' stellar mass and velocity dispersion. iMagorrian et al .111 9981 : iFerrarese fc Merrittil2000t iGebhardt et al.ll2000 ) 



imply that this connection spans an astonishing range of spatial scales; the AU-sized accretion disk whose 
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characteristics determine the black hole growth rate must receive information about the conditions of the 
galaxy on scales of kiloparsecs, or eight orders of magnitude larger. These correlations between black hole 
mass and various bulge properties hint at simultaneous evolution of AGNs and their host galaxies. 

The processes that are involved in the maintenance of the M — a relation may also govern the evolution 
of host galaxy morphology and color. In the local universe, galaxies exhibit bimodal distribu tions of several 



critical parameters, including rest-frame optical color, star formation rate, and gas fraction (jStrateva et al 



20011 : iKauffmann et al.l l2003t iBaldrv et al.l |2004 etc) . Star formation must be suppressed significantly to 



explain the low star formation rates and the redness of early- type, passively evolving galaxies today. There 
is circumstantial evidence for the participation of AGN in this suppression of star formation, or "quench- 
ing," including their presence in the green vall ey of the galaxy color-magnitude diagram (CMP), intermedi 



ate between the blue cloud a nd red sequence (JKauffmann et al 



Xue et al 



201 



w. 



2003; 



Nandra et al. 



2007 : 



Schawinski et al 



2009 : Cardamone et al.ll2010l ). although the AGN fraction does not appear to be higher in the green valley 



Major mergers of gas-rich galaxies can result in the rapid growth of central black holes. In this pic- 
ture, the co llision of gas clo uds saps angular momentum, efficiently sending gas to the central regions of 
the galaxy fflernq uisti Il989 ). where it is largely consumed in star formation and partially accreted onto 
one or more BHs. In several models, the "M-cr" relation is m aintained through self-regulating processes 
involved in galaxy a nd supermassive black hole (SMBH) growth ( Hopkins et al.ll2005t iDi Matteo et al.ll2005t 
Hopkins et al.ll2009l and references within). This self-regulation may be modulated by star formation, in 
which mo mentum iniection by Type la supernovae and radiation pressure from O stars disrupt further star 
formation (jMurrav et al.ll2005[ ). However, the immense energies associated with gravitational accretion onto 
a SMBH dominate the energy budget of the galaxy; this suggests that the AGN itself may provide sufficient 



energy to remove gas and limit the growth of the BH (ISilk fc E 



lees 



1998 



Hachnelt et al.lll998t fcnsjbooj 



Pi Matteo et allboosl : iHopkins et alTboOst ISpringel et al.lbood : iHopkins et al...2006: iMenci et al.ll2008^ ~A 



further conseque nce of the most energetic AGN feedback processes may be the e xpulsion of a majority of 
the galactic gas (jSanders et al.lll988t |Pi Matteo et al.ll2005t IHopkins et al.l 120061 1. This quenching of star 
formation may be a mechanism by which AGN activity controls the evolution of galaxies from the blue cloud 
to the red sequence. 

One feature of a model that includes rapid, AGN- induced quenching is an evolution between obscured 
and unobscured AGN states (JHopkins et al.ll2006t iMenci et al.l l2008) . The ejection of gas/dust in the central 
core of the galaxy is quickly followed by gas removal on kpc-scales. This process naturally produces a positive 
correlation between AGN ob scuration and the mass fraction of young stars (a consequence of sufficient gas 
supply) in the host galaxies (JMenci et al.l 120081 ) . Unobscured AGN are expected to reside in relatively gas- 
free host galaxies, cleared by a previous blowout phase. These hosts are expected to be redder than the 
hosts of obscured AGN, which would have higher star-formation rates. 

It is not clear how this theoretical evolution between obscured and unobscured states is consistent with 
the observed AGN unification paradigm. AGN unification schemes attempt to explain the properties of the 
many observed types of AGN with a common physical model. The most successful unification picture pos- 
tulates that the differences between Seyfert types in t he local universe are due to variati on of the orientation 
of the observer relative to a centralized dusty torus (|Antonuccil Il993l : lUrrv et al.l 119951) . Since obscuration 
is solely determined by geometry in a unified scheme, there is no reason to expect differences between the 
properties of the host galaxies of obscured and unobscured AGN. 



Many comparisons of obscured and unobscured AGNs have been performed at low and high redshift, 



using both optical- and X-ray-based diagnostics of AGN obscuration and various measures of star formation 
rate or gas content. Locally, obscured Seyferts a re associat ed with sign ificant star f orming activity, in some 
cases mo re than compa rable unobscured hosts ( Cid Fernandes fc Terlevichlll995t iGonzalez Delgado et al. 



20011 ) ■ iLacv et al.l (|2007l ) finds that obscured Type II QSOs at redshifts of 0.3 < z < 0.8 are associated with 
significant star forming activity. Howe ver, using HI measurements, unobscure d broad-line AGN hosts at 
low-redshift are observed to be gas- rich ( Ho et al.ll2008 ). ISchawinski et al.l ( 20091) find no dependence of host 



galax y color on obscuring col umn for a small sa mple of SWIFT BAT AGN (JGehrels et al.ll2004i iTueller et al 



20081 ) in the local universe. iPage et al.l (|2004l ) find stronger submillimeter emission in an X-ray obscured 



sample compared t o an X -ray unobscured sample, implying a higher star formation rate in the obscured 
sample. iGilli et al.l ( 2009fl find that the clustering lengths of z ^ 1 X-ray selected AGNs with broad-lines 
and those that lack broad lines are similar. They also note that the clustering length of the X-ray hard 
sources is statistically simila r to that of the soft sources, although the errors are dominated by number 



Shao et al.l (|2010l ) measure the star formation rates of X-ray selected AGNs in GOODS-North 



statistics. 

with Herschel far-infrared fluxes, finding no correlation with X-ray absorbing column. 



Pierce et al.l (|2010bf ). GPIO hereafter, uses optical colors of X-ray selected AGNs at z ~ 0.7 in AEGIS 
to compare hard and soft sources. Using colors measured with both integrated and extended apertures, they 
find that unobscured sources are bluer than obscured sources in rest-frame U — V. 



1.1. Comparing the Colors and Color Gradients of Obscured and Unobscured AGN 

Generally, understanding the relationship between AGN host galaxy properties on kpc-scales and the 
properties of the AGN on smaller scales can be helpful in constraining the role of AGN in the evolution 
of galaxies. We focus on comparing the host properties of obscured and unobscured sources, given the 
hypothesis that unobscured sources will reside in redder hosts in a feedback scenario involving AGN-driven 
blowouts. We adopt a two-pronged approach: We utilize two imaging datasets of different spatial resolution, 
allowing us to probe host properties on two spatial scales. We first measure the rest-frame u — R colors of 
a large sample of AGN host galaxies at high galactic radii (r ^ 6 — 12 kpc) and search for correlations of 
color with AGN obscuration or hardness ra tio. The combination of opti cal and infrared colors is a powerful 
diagnostic of the presence of young stars ( Gil de Paz fc Madord 120021) . particularly rest-frame u — R, or 
observed V — J at z ^ 0.8. This sample is limited in spatial resolution by the infrared imaging data (O'.'S - 
0'.'6 for ISAAC imaging and - 0"25 for NIC3). 

We then turn to a smaller sample of AGN host galaxies imaged by the Wide Field Camera 3 In- 
frared Channel (WFC3 /IR) on HST as part of the Early Release Science (ERS) demonstration program 
( Windhorst et al.ll2010l) . The improved spatial resolution of ~ 0'.'2 enabled by WFC3 allows us to probe 
rest-frame NUV — R color closer to the central regions of host galaxies at r ^ 1 kpc, where AGN-induced 
winds may be the most energetic and color effects potentially more drastic. 

Section 2 identifies th e VLT /ISAAC and AC S data we use in the Great Observatories Origins Deep Sur- 
vey (GOODS) field South (^ Giavalisco et al.ll2004 ) and the sample selection. Section 3 details the methodology 
used to obtain uncontaminated aperture photometry of AGN host galaxies in these fields. Section 4 describes 
a suite of cloning simulations used to check the analysis and derive error bars. Section 5 presents results 
for all samples and Section 6 concludes. We assume a flat cosmology throughout, with Hq = 70 km/s/Mpc, 
riM = 0.3, and 17a == 0.7. 



DATA 



This paper focuses on a set of ^ 78 X-ray s elected sourc es in the Great Observa tories Origins Deep Survey 



(GOODS) field South (iGiavalisco et alj 120041 ) and North (1 Alexander et alj 120031 ). overlapping the Chandra 
Deep Fields (CDF). The deep, 1 Megasecond Chandra exposure i n the CDF-S and the 2 Ms exposure in 
the CDF-North hav e revealed hundreds of AGNs at0.5<z<1.5 |Giacconi et al.ll2002l : Barger et al.ll2003 : 
Szokolv et al.l 120041 ). Deep optical and near-infrared imaging from the HST Advanced Camera for Surveys 
(ACS) is available in four filters, F435W (B), F606W (V), F775W (i), and F850LP (z). This paper makes 
use of deep J,H, and Ks imaging in t he CDFS usin g the ISAAC instrument at the VLT. The data reduction 
and depths are described in RctzlafF et al.l (|2010bf ) . 



We utiliz e deep FllOW and F1 60W imaging publicly available in the GOODS-North and South fields. 
These NIC-3 ([Thompson et al.lll998l ) pointings have been acquired through a variety of programs, including 
Guest Observer, DDT and Parallel ob serving programs, to taling 1400 orbits of FllOM^ and F160W imaging. 
As part of an HST archival program, iMagee et al.l (|2007l ) ha ve reduced and mos aiced all NIC-3 imaging in 
GOODS. This data reduction procedure is described fully in Magee et al.| (12003). Further details regarding 
the use of this NIC-3 data set are given in section 4.2.1 of lAmmons et al.l (|2009l ). 

The ne ar-ultraviolet (NUV) filter is defined a s used for the GALE X mission ( Martin et al.l2005tlMorrissev et al 
20051120071) . The NUV - R color is defined as in lSalim et al.l (J2007t ). All magnitudes are in AB units unless 



explicitly stated otherwise. Selection criteria are enumerated in section [23 



2.1. Near-Infrared ISAAC Imaging 

The ISAAC instrument on the VLT has been used to mosaic the ACS region of GOODS-So uth to median 
5(T point source depths of Jab = 25.2, Hab = 24.7, and Ksab = 24.7 JRetzlaff et al.ll2010bl) . The spatial 
coverage in each of the bands is 172 A, 159.6, and 173.1 arcmin^ in J, H, and Ks respectively. The exposure 
times varied from 3-8 hours in each band, totalling 1.3 Megaseconds. This program was queue-scheduled 
to optimize image quality; the image FWHM's range from 0'.'3 to 0^.'65. Fully reduce d and photometrically 
calibrated data is available from the GOODS ISAAC webpage (JRetzlaff et al.ll2010a[ ). 



2.2. Near-Infrared NICMOS Imaging 

In this study, we use nearly 100 square arcminutes of NIC-3 FllOW and F1601^ imaging in GOODS 
North and South. These pointings, which have been acquired thro ugh a variety of ob serving programs 
and are publicly available, have been fully reduced and mosaiced by iMagee et al.l (J2007l ). This reduction 
procedure includes basic calibration (removing the instrumental signature), cosmic-ray removal, treatment 
for post-SAA cosmic ray persistence and electronic ghosts, sky subtraction, non-linear count-rate correction, 
artifact masking, robust a lignment and regist ration for large mosaics, weight map generation, and drizzling 
onto a final image frame (JMagee et al.l 120071 ). Sensitivities vary across the GOODS fields due to differing 
levels of exposure. 



2.3. WrC3/IR ERS Imaging 

We use public WFC3/IR imaging in the ERS region with the F125W and F160W filters. The data 
were reduce d using the standard c alwfS pipeline. After pipeline calibration, a residual gain "pedestal" was 
noticed f see IWindhorst et al.ll201 CF) and corrected through a multiplicative scaling of each quadrant. Images 
were corrected for geometric distortion and drizzled onto a 90mas grid using the MultiDrizzle package 



Berti 



Koekoemer et al.ll2002L 120071 ). Individual pointings were registered to a single astrometric grid using SWarp 



20021 ). which corrected for a small rotation offset in the astrometric solution from MultiDrizzle. 
A final mosaic of the data was created with SWarp and a small residual sky gradient was removed. The 
compact source sensitivity is 27.5 and 27.2 magnitudes AB for a 5ct detection in the F125W and F160H^ 
bands, respectively. 



2.4. Sample Selection 

The overall sample is composed of X-ray selected sources in the GOODS ACS regions with redshifts, as 
described below. The overall sample is divided into subsamples for further analysis with cuts in the X-ray 
hardness ratio and X-ray luminosity. The samples and selection criteria are enumerated in Table [TJ 

This sample i s X-ra y selected using the 1 Ms Chandr a Deep Field South (C DFS) exposure as described 



in lGiacconiet al.l (|2002[ ) and the 2 Ms CDF-N exposure (lAlexander et al.ll2003D . Ah GOODS-South X-ray 
sources were chosen with a SExtractor ( Bertin fc ArnoutslllQQa) S/N detection thre shold of 2.1. GOODS- 



North X-ray point sources are identified with WAVDETECT ( Freeman et al.l 120021 1 using a false-positive 
probability threshold of 10^^ in one of seven X-ray bands ranging from 0.5 keV to 8.0 keV. This threshold 
is comparable to the GOODS-South S/N detection threshold of 2.1. 

Of 339 X-ray sources detected in GOODS-South i Giacconi et al.l 120021 ) . 191 unique sources are identi- 
fied as extragalact ic (non-stellar), he in the ACS and IS AAC tiled region, and have photometric redshifts 



from COMBO- 17 (I Wolf et al. 



2004 



Zheng etal.1 120041 ) or iLuo et al.l (|2010l ) or spectroscopic redshifts from 



a 



Szokolv et al.l ( 20041 ) or lSilverman et al.l ( 2010l ). Sources with extended X-ray emission are omitted from the 
sample. A redshift cut isolating sources with 0.5 < z < 1.5 and an optical magnitude cut of i? < 24 trim the 
sample to 56 objects. 



Of the 503 X-ray p oint sources d e tected in GOODS-North ([Alexander et al.l l2003l ) . 284 have reliable 
spectroscopic redshifts in lBarger et al.l ( 20031) . 22 sources possess spectroscopic redshifts, fall in the redshift 
range 0.5 < z < 1.5, appear within the NIC-3 F160W tiled region in GOODS-North, have Lx,o.5-skeV > 10"^ 
ergs s^^, and have R < 24. The median X-ray detection S/N between 0.5 and 8.0 keV for the remaining 
sample is 11.95 and the minimum is 4.91, above the detection threshold for the GOODS-South catalog. 
Including AGN from both northern and southern fields, 78 AGN remain in the sample. The identification 
numbers, redshifts. X-ray luminosities, hardness ratios, u — R colors, and U ~ B colors are given in Tabled 

92% of the sample of 78 AGN hosts have secure spectroscopic redshifts. We divide the sample at a 
hardness ratio of —0.2 to isolate obscured and unobscured sources. Sources with HR < —0.2 are termed 
unobscured or soft and sources with HR > —0.2 are termed obscured or hard. These labels are used to refer to 
th ese types of sou r ces th roug hout this paper We u se the hardness ratio as defined as HR — {H~S)/{H + R) 



Giacconi et al.l (|2002l ) and ISzokolv et al.l (|2004l ). where H and S are the net instrument count rates in the 



hard (2-10 keV) and the soft (0.5-2 keV) bands, respectively. 

Hardness ratio is expected to be an imperfect measure of AGN obscuration due to K-correction effects. 



We assess this usi ng a subsample of GOODS-South X-ray sources with full X-ray spectral fits available in 
Tozzi et al.l (|2006l ). A hardness ratio of —0.2 at z ~ 1 corresponds to a neutral hydrogen column density of 
2 X 10^^ cni~^. The neutral hydrogen colu mn corresponding to HR = —0.2 shifts by a factor of ^ 2.5 from 
z ^ 0.5 to z '-^ 1.5 due to the k-correction (JTozzi et al.ll2006l ). The range of hardness ratios that correspond 



to 2 X 10^^ cm ^ over the redshift range 0.5 < z < 1.5 is —0.3 < HR < —0.1. 25% of soft sources selected 
with HR < -0.2 have Na > 2x 10^^ cm'^; 18% of hard sources selected with HR > -0.2 have Nh < 2x 10^^ 
cm~^. This low level of sample mixing does not affect the conclusions of this paper. 

The final redshift ranges, mean redshifts, and numbers of sources in each subsample are indicated in 
Table [TJ Table [1] lists three groups, each of which are further subdivided into hard and soft subsamples: 
"all", "lum", and "ERS." The "all" sample is selected according to the criteria listed above. The "lum" 
group is a subsample of the "all" group and only includes objects more luminous than Lxfi.s-iokeV — lO"'^'^ 
ergs s~^. The "ERS" subsample only includes sources in the ERS region sampled by HST WFC3 imaging. 

To enable unbiased comparison of obscured and unobscured host properties, we cut the redshift ranges 
of the hard and soft subsamples to impose the mean redshifts of these two samples to match. The redshift 
distribution of soft sources is generally weighted towards higher redshifts, so we trim the highest redshift 
sources from the soft subsamples. Similarly, we trim low redshift sources from the hard samples. As a result, 
both the mean redshifts and the mean X-ray luminosities of the soft and hard subsamples within each group 
match closely. As shown in Figure [U the trimmed redshift histograms of the two populations are broadly 
similar. Note in Table [T] that both the mean redshift and mean X-ray luminosity of the "lum" group is 
larger than the "all" group, as expected for a subsample of more luminous objects. The mean redshift of 
the "ERS" group is significantly lower than that of the "all" group because the ERS region overlaps with a 
known galaxy overdensity at z ^ 0.73. The mean X-ray luminosity of the "ERS-hard" subsample and the 
"ERS-soft" subsample are significantly different, 10^'^°^ ergs s^^ and lO^'^-^^ ergs s~^, respectively; this is 
due to two particularly X-ray luminous sources in the small "ERS-soft" subsample. If these are removed, 
the mean colors of the "ERS-soft" sample do not change more than the la measurement error. 



Table 2:: Tabulated redshifts. X-ray luminosities, hardness ra- 
tios, and measured colors for individual sources. Sources identified 



with "LUO" are GOODS-South sources as numbered in Luoetal 



(|201flf ). Sources identified with "XID " are GOODS-Sou th sources 



as num bered in ISzokolv et al.l (|2004( ) but not found in iLuo et al. 



(|2010l ). Sour ces identified with "ABB " are GOODS-North sources 
num bered in [Alexander et al.l (|2003l) with spectroscopic redshifts 
from lBarger et al.l (J2003l ). 



XID 


z 


type 


Lx 


HR 


yU f^) outer 


{U — Bjouter 


{u - r)int 


(U - B)„t 




LUO 408 


1.23 




43.79 


-0.45 


1.66 


0.564 


1.70 


0.564 


LUO 370 


1.02 




43.48 


0.12 


1.87 


1.04 


1.71 


0.694 


LUO 335 


1.22 




43.46 


-0.2 


2.18 


0.695 


0.341 


0.0121 


LUO 324 


0.84 




44.03 


-0.53 


2.07 


0.628 


0.683 


0.0259 


LUO 319 


0.66 




43.04 


-0.49 


2.48 


0.716 


1.39 


0.435 


LUO 316 


0.67 




43.36 


-0.4 


1.93 


0.723 


1.99 


0.665 
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XID 


z 


type 


Lx 


HR 


(m - r)outer 


(U — B)outer 


{u - r)^nt 


{U - BUt 




LUO 302 


0.84 


1 


43.21 


-0.32 


1.09 


0.559 


1.46 


0.529 


LUO 288 


1.03 


0.5 


43.44 


-0.37 


1.48 


0.521 


1.78 


0.582 


LUO 283 


0.96 


0.5 


43 


-0.23 


0.944 


0.225 


1.28 


0.341 


LUO 273 


0.74 




43.52 


-0.56 


1.64 


0.679 


0.722 


0.241 


LUO 267 


1.22 




44.19 


-0.47 


1.83 


0.669 


0.459 


0.0749 


LUO 249 


0.67 




43.29 


0.52 


1.77 


0.566 


1.82 


0.669 


LUO 247 


0.73 




44.48 


-0.54 


1.83 


0.669 


0.380 


0.0914 


LUO 246 


0.73 




43 


0.06 


2.00 


0.669 


2.01 


0.662 


LUO 242 


1.03 




44.23 


-0.63 


1.08 


0.454 


0.386 


0.135 


LUO 224 


0.73 




43.23 


0.44 


2.07 


0.771 


1.67 


0.609 


LUO 216 


0.53 




42.78 


-0.47 


2.23 


0.657 


2.38 


0.704 


LUO 168 


1.1 




44.21 


0.61 


1.65 


0.755 


1.91 


0.698 


LUO 167 


0.57 




43.17 


-0.55 


1.91 


0.588 


1.74 


0.546 


LUO 145 


0.67 




42.92 


-0.45 


2.24 


0.794 


2.20 


0.783 


LUO 127 


0.61 




43.42 


0.11 


1.88 


0.589 


1.45 


0.470 


LUO 49 


1.04 




43.73 


-0.44 


2.22 


0.789 


0.620 


0.168 


LUO 270 


0.96 




43.32 


-0.54 


1.43 


0.644 


1.45 


0.458 


LUO 311 


1.31 




42.64 


-1 


0.316 


0.0924 


0.668 


0.168 


LUO 118 


1.03 




42.8 


0.1 


1.89 


0.715 


1.74 


0.613 


LUO 228 


1.09 




43.56 


1 


1.90 


0.764 


1.99 


0.752 


LUO 189 


0.6 




43.3 


1 


1.86 


0.643 


2.28 


0.679 


LUO 90 


0.55 




42.49 


0.16 


1.61 


0.604 


1.78 


0.583 


LUO 381 


0.66 




42.32 


-0.17 


2.07 


0.719 


2.20 


0.752 


LUO 304 


0.84 




42.84 


0.01 


1.13 


0.271 


1.57 


0.333 


LUO 415 


1.14 




42.42 


-1 


1.92 


0.712 


2.11 


0.757 


LUO 393 


0.67 




43.9 


0.39 


2.39 


0.833 


2.27 


0.774 


LUO 109 


0.93 


0.9 


42.83 


0.13 


1.35 


0.531 


1.67 


0.619 


LUO 204 


0.73 




42.38 


1 


2.09 


0.766 


2.46 


0.796 


LUO 310 


0.73 




43.22 


1 


2.09 


0.748 


2.05 


0.717 


LUO 211 


1.22 




42.65 


-1 


1.08 


0.374 


1.45 


0.462 


LUO 254 


0.74 




42.19 


-1 


1.99 


0.654 


2.29 


0.804 


LUO 115 


0.76 


0.6 


42.17 


-0.47 


1.44 


0.331 


2.03 


0.451 


LUO 386 


1.17 




43.43 


0.52 


1.11 


0.462 


1.71 


0.579 


LUO 226 


1.04 




43.23 


0.23 


1.52 


0.746 


1.27 


0.554 


LUO 260 


1.32 




43.45 


0.6 


1.08 


0.289 


0.834 


0.244 
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XID 


z 


type 


Lx 


HR 


(m - r)outer 


(U — B)outer 


{u - r)int 


{U - BUt 




LUO 131 


0.74 




43.54 


1 


1.61 


0.599 


1.89 


0.631 


LUO 72 


0.72 




43.39 


1 


1.78 


0.377 


1.75 


0.357 


XID 511 


0.77 




42.01 


-0.24 


0.868 


0.379 


0.952 


0.411 


LUO 298 


0.67 




42 


-0.14 


1.98 


0.694 


2.13 


0.743 


XID 516 


0.67 




42.36 


-0.26 


0.726 


0.330 


1.02 


0.371 


LUO 235 


1.03 




42.62 


-0.03 


0.930 


0.452 


1.05 


0.419 


LUO 178 


0.96 




42.46 


0.18 


1.93 


0.667 


1.85 


0.642 


LUO 117 


0.68 




42.49 


0.34 


2.14 


0.748 


2.17 


0.774 


LUO 114 


0.58 




42.41 


-0.01 


1.72 


0.498 


2.00 


0.632 


XID 580 


0.66 




42.32 


-1 


1.90 


0.730 


2.21 


0.726 


LUO 129 


1.33 




43.72 




1.61 


0.472 


1.95 


0.551 


LUO 177 


0.74 




44.37 




2.03 


0.821 


2.28 


0.839 


LUO 195 


0.67 




43.3 




2.01 


0.694 


2.28 


0.777 


XID 611 


0.98 




43.41 




1.38 


0.646 


2.00 


0.635 


XID 612 


0.74 




43.3 




1.49 


0.527 


1.68 


0.570 


ABB 40 


1.379 




44.42 


-0.43 


1.42 


0.299 


1.24 


0.228 


ABB 55 


0.64 




42.04 


-0.1 


1.67 


0.567 


1.71 


0.576 


ABB 82 


0.68 




42.6 


0.34 


1.73 


0.642 


1.93 


0.723 


ABB 90 


1.14 




43.42 


0.5 


1.57 


0.575 


1.72 


0.504 


ABB 113 


0.85 




43.2 


-0.39 


2.01 


0.783 


2.19 


0.670 


ABB 115 


0.68 




43.49 


-0.34 


1.71 


0.663 


1.19 


0.423 


ABB 121 


0.52 




42.35 


0.2 


1.73 


0.476 


1.92 


0.524 


ABB 142 


0.75 




42.61 


-0.07 


2.21 


0.814 


2.22 


0.775 


ABB 150 


0.76 




42.89 


0.49 


1.50 


0.582 


1.79 


0.667 


ABB 157 


1.26 




43.99 


0.38 


1.15 


0.341 


1.31 


0.365 


ABB 158 


1.01 




43.11 


0.34 


2.09 


0.725 


2.23 


0.737 


ABB 177 


1.02 




42.66 


-0.32 


1.19 


0.432 


1.54 


0.504 


ABB 193 


0.96 




43.67 


-0.55 


0.878 


0.362 


0.742 


0.279 


ABB 201 


1.02 




42.8 


0.64 


1.91 


0.786 


1.73 


0.601 


ABB 205 


1.38 


0.5 


43.51 


-0.33 


1.21 


0.259 


1.55 


0.314 


ABB 212 


0.94 




42.26 


-0.37 


2.19 


1.08 


1.80 


0.756 


ABB 217 


0.52 




42.2 


0.68 


2.63 


0.768 


2.51 


0.721 


ABB 222 


0.86 




42.93 


-0.22 


1.90 


0.654 


1.96 


0.664 


ABB 351 


0.94 




42.14 


-0.3 


1.77 


0.650 


1.99 


0.680 


ABB 352 


0.94 




42.55 


0.53 


1.49 


0.547 


1.61 


0.587 
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XID 


z 


type 


Lx 


HR 


[u - r)outer 


{U — B)outer 


{u - r),nt 


{U - BUt 




ABB 384 


1.02 


1 


43.52 


0.77 


1.74 


0.648 


1.95 


0.675 


ABB 451 


0.84 


1 


44.04 


-0.42 


1.44 


0.568 


0.893 


0.260 





2.4.1. Color Magnitude Diagram 

Figure [2] plots the rest-frame u~ R vs Mu color magnitude diagrams (CMD) of the "all" sample. The 
rest-frame photometry is obtained by performing 3'.' filled-aperture photometry on all optical/NIR bands for 
all AGN, interpolating linearly through the SED, and sampling the de-redshifted SED at the appropriate filter 
wavelengths. For comparison, rest-frame u — R and Mr measurements for the underlying galaxy population 
are plotted as small black circles for qualitative compar ison. This sample of 535 control galaxies is selected 



from the GOODS MUSIC catalog (jSantini et al.l 120091 ) identically to the AGN sample (i.e., spectroscopic 
redshifts with 0.5 < z < 1.5 and R < 24) except that X-ray sources with Lx > 10"*^ ergs s~^ are excluded. 
Rest-frame colors are computed via interpolation through observed ACS/ISAAC bvizJHK photometry, as 
described in Section [31 A ^ 0.15 magnitude correction is added to the IR bands to correct for a difference 
in photometric aperture size between this sample and the MUSIC sample. 

In Figure [21 it is apparent that the AGN tend to reside in the most luminous systems . Both red 



seque nce galaxies (u — i? ~ 2.0) and blue cloud galaxies {u — R^ 1.0) are present in this sample (jSalim et al. 



20071 ). The mean loci of the soft and hard sources are shown as blue and red error bars, respectively. 
Using integrated photometry, it immediately appears that the soft sources are significantly bluer and more 
optically luminous than the hard sources. However, integrated photometry is likely contaminated by the 
nonthermal emission produced by the AGN engine. The next section introduces our method of measuring 
uncontaminated rest-frame host colors. 



3. APERTURE PHOTOMETRY 

We now measure aperture photometry with an elliptical annular aperture with inner de-projected radius 
of 0'.'75 and outer de-projected radius of 1'.'5, hereafter labeled the "outer" aperture. We select a large inner 
radius of 0'.'75 that minimizes contamination from central point sources. This dimension is oversized with 
respect to the spatial resolution of the near-IR imaging, which is set by ground-based seeing. The outer 
radius of 1'.'5 is chosen to maximize the total signal-to-noise ratio over all bands for typical galaxy profiles 
of sersic index n = 2. 

The axis ratios and position angles of the annuli are determined for each AGN host separately by 
manually fitting to isophotes in the ACS z-band image. This is performed by lining up ellipses with major 
axis lengths of 3'.'0 with variable axis ratios and position angles until the nearest isophote is visually matched. 
Only one annular aperture is used for each AGN host. 

Rest-frame colors are computed from observed photometry with spline interpolations using IDL's interpol 
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Fig. 1. — Redshift histograms for soft (blue line) and hard (red line) sources. Soft and hard sources are 
divided at a hardness ratio of —0.2. The mean redshifts of the two populations are shown as dashed lines. 
The mean redshifts of the samples match when the redshift limits are trimmed as described in Section 12.41 
The redshift distributions of the soft and hard sources are broadly similar, both displaying peaks at known 
redshifts of 2; ^ 0.7 and z ^ 1.0. 
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Table 1: Table of sample sizes, redshift ranges, mean redshift, mean X-ray luminosities, summarized color 
information, and selection criteria for all subsamples. 



Sample 


# 


z range 


(^> 


(i^) 


((u-r)i„t) 


((W - r)outer) 


Selection 




all 


78 


0.5-1.5 


0.877 


43.11 


1.65±0.09 


1.69±0.10 


Lx>42 




all-hard 


42 


0.55-1.5 


0.869 


43.10 


1.81±0.10 


1.74±0.11 


Lx > 42, HR>- 


-0.2 


all-soft 


31 


0.5-1.3 


0.864 


43.12 


1.45±0.11 


1.65±0.12 


Lx > 42, HR<- 


-0.2 


lum-all 


19 


0.5-1.5 


1.018 


43.92 


1.31±0.12 


1.65±0.14 


Lx > 43.5 




lum-hard 


8 


0.5-1.5 


0.994 


43.85 


1.94±0.17 


1.76±0.18 


Lx > 43.5, HR > 


-0.2 


lum-soft 


10 


0.5-1.3795 


1.001 


44.01 


0.78±0.16 


1.61±0.17 


Lx > 43.5, HR < 


-0.2 


ERS-all 


24 


0.5-1.5 


0.788 


43.21 


1.61±0.11 


1.78±0.13 


Lx>42 




ERS-hard 


15 


0.5-1.5 


0.793 


43.02 


1.85±0.15 


1.74±0.16 


Lx > 42, HR>- 


-0.2 


ERS-soft 


9 


0.5-1.5 


0.781 


43.54 


1.20±0.16 


1.85±0.17 


Lx > 42, HR<- 


-0.2 
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Fig. 2. — Rest frame AB u ~ R vs. Mn color-magnitude diagram of the "all" sample, using integrated 
photometry. Large circles denote strong {Lxfi.b-iokeV > 10'*'^ ergs s"^) sources and small circles denote weak 
{Lx,o.5-iokeV < 10'*'^ ergs s"-^) sources. Hard, obscured sources (HR > —0.2) are shown with overplotted 
squares and red colors. Soft sources (HR < —0.2) are shown with blue colors. Rest-frame values are 
estimated from observed photometry as in the text. Both u — R colors and M^ are measured using a 
filled 3'.' diameter circular aperture. Colors and absolute magnitudes are uncorrected for central point source 
contamination. Blue and red error bars denote the mean colors and luminosities for the soft and hard sources, 
respectively, with errors calculated as in Section 15.41 Small black circles denote the colors and magnitudes 
of the underlying galaxy population, selected as explained in the text. 
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function and known filter transmission functions. The spline-fit SED is constrained to match the observed 
photometry upon integration with the appropriate filter curves. Rest-frame magnitudes in u, R, and other 
bands are computed by integrating the de-redshifted SED with these filter functions. We avoid using galaxy 
template spectra for fitting photometry, which can be problematic with the addition of non-thermal AGN 
continuum emission. We test our techniq ue by compar i ng res t-frame u — R colors measured in a filled 2'.'5 
diameter aperture with those computed bv lTavlor et al.l ( 20091) using SED-fitting techniques. For 44 isolated 



sources in GOODS-South that overlap the two samples, the u — R colors match with a symmetric scatter 
of cr = 0.15 magnitudes and a minor systematic offset of 0.015 magnitudes. The scatter shrinks to cr = 0.11 
magnitudes for sources not obviously dominated by central point sources, similar to the average level of 
measurement error in the sample. We conclude that our technique is appropriate for measuring the average 
rest-frame colors of subsamples of galaxies. 

The "outer" aperture is expected to be robust to contamination from central point sources associated 
with nonstellar AGN continua. However, light from the central AGN may be scattered into the "outer" 
aperture by telescope diffraction, or the effect of a non-ideal Point Spread Function (PSF). PSF contamination 
can be non-trivial for cases in which the AGN emission is dominant over the light from the host galaxy. 
These cases require either a correction to the aperture photometry or special image processing to reveal 
the host galaxy. We present two techniques for correcting aperture photometry for the effect of scattered 
AGN contamination. Firstly, Section IXTj describes a method of estimating contamination from the PSF and 
subtracting its contribution from measured phot ometry. Secondly, Section [3.21 describes the implementation 



of the CLEAN (|Hogbomlll974t iKeel et al.lllQQll ) deconvolution algorithm to remove contamination directly 



from the images before measuring photometry. Both methods requires excellent constraint of the PSF. 



3.1. Correcting "Outer" Aperture Colors for PSF Contamination 

Central embedded point sources contribute flux to annular photometric apertures via PSF scattering. 
We approximate the contaminating flux Fcontam for a given source in a single band S as 

/ , FpsF,S 

P _ outer V^ p 

^ contam — ^ — ^ ^ / ^ ^ imagers 

/ , FpsF.S core 

core 

where Fimage,s is the flux distribution of the source in band S and FpsF.s is the flux distribution of 
the centered PSF in band S. Y,outer denotes a sum of fluxes over the "outer" annular aperture. 'S'core here 
denotes a sum over a filled, circular aperture with radius O'.'IS. Note that the first fraction is similar to a 
concentration measurement, except that wider aperture is elliptical and annular. The decontaminated flux 
is then given by 



^ decontam I / ^ ^ imagers I ^ contani 

\outer / 

This method corrects for all PSF-scattered flux originating from within the central core aperture, no 
matter the source. PSF-scattered flux originating from outside of the core aperture is not removed, but it is 
assumed that this light is intrinsically stellar. 
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3.2. CLEAN Method 



This section describes the use of the CLEAN algorithm to deconvolve galaxy images and allow estimation 
of intrinsic, uncontaminated photometry in a variety of photometric apertures. 

The CLEAN algorithm was originally designed to rem ove confusing sidelobe patterns in synthesized- 
beam radio interf erometry at high resolution ( Hogbomlll974 ) but has been modified to perform deconvoluti 



ion 



of direct imaging ( Keel et al.lll991l ). CLEAN iteratively constructs a model of the instrinsic light distribution, 
returning both the model and the image residual. In each iteration, CLEAN adds "points" (delta functions) 
to the model that correspond to the location of the maximum residual in the original image. The intensity 
of the point is a fraction (gain) of the intensity of the maximum residual. The convolution of the PSF with 
that point is subtracted from the original image. The loop repeats until all of the flux from a specifled region 
in the measured image has been remove d or after a max imum number of iterations. CLEAN preserves the 



photometric scale of the original image (|Keel et al.lll99l[ ). 



For this study, CLEAN is performed over a 4'.'5 box in both ACS and ISAAC images. The chosen 
number of maximum iterations is 500 and the gain is 5%. The algorithm is modified to find peaks using a 
cross-correlation of the PSF against the image rather than subtracting at the maximum pixel location. 

The deconvolved image is the sum of the residual image and the CLEAN model. Galaxy fiuxes uncon- 
taminated by the wings of any central point source can then be measured by performing aperture photometry 
directly on the deconvolved image with a specified aperture. 



3.3. Stacking Analyses 

Stacking sky-subtracted, reduced images of galaxies in various bands is helpful for revealing the mean 
flirxes and colors of various sub-populations, and serving as a check on other methods. In this study, we 
stack images that have been masked to remove contaminating background sources and deconvolved with 
CLEAN. We use observed V and J colors, which correspond to nearly u and R at z ^ 0.8. The gain in S/N 
realized by stacking is tempered by two major disadvantages: (1) Observed-frame colors do not represent 
the mean rest-frame colors, as the sources lie at a variety of redshifts; and (2) Stacking requires the use of 
a single photometric aperture on a stacked image rather than a variety of apertures that correspond to the 
position angle and axis ratio of individual sources. With these caveats in mind, we use a circular, annular 
photometric aperture of inner radius 0'.'75 and outer radius I'.'S to compute observed V — J colors. 



3.4. Sky Subtraction 

Meaurements of surface brightness at high galactic radii are highly sensitive to errors in sky subtraction. 
The sky value for each band is measured by manually masking 3-sigma sources and computing the geometric 
mean of the values in a circular annulus with inner radius 4'.'0 and outer radius 6'.'0. These radii were 
determined empirically by selecting empty regions in ACS and ISAAC fields and attempting to predict the 
sky value in an inner circle of S'.'O diameter; this combination of sky radii minimized the errors in all ACS 
and ISAAC bands. It is important to fix the annular radii across bands so as not to measure artificial colors 
from galactic contamination. 

Other methods of measuring sky were tested, including using sky measurements at varying radii to 
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extrapolate the sky value inwards to zero radius. The simple masked mean is superior to extrapolation. 



3.5. PSF Selection 

PSF reference stars are used to correct the annular photometry for the presence of bright, centralized 
point sources. It is important that these PSF references be reliable to enable good recovery of underlying 
galaxy photometry. 

The ACS, NIC3, and WFC3 point spread functions are spatially and temporally stable compared to 
ground-based PSFs. For each ACS, NIC3, or WFC3 band, images of bright unsaturated stars are taken 
from locations in the field. The radial profiles of the stars are compared to Tiny Tim PSFs to verify that the 
sources are unresolved. The PSFs are chosen to be in empty, isolated regions; any sources in the area are 
masked with replicated boxes of pixels taken from nearby empty regions with the area. A sky subtraction 
is performed using the masked mean method described in the above section. The PSFs are centered by 
sub-pixel interpolation. 

For ISAAC PSF reference stars, the sources are verified as unresolved by comparing to ACS imaging. In 
each ISAAC field, the brightest, isolated star with pixel values less than 10 times saturation level is selected. 
Sky subtraction and centering are performed on these PSFs as described above. The ACS and WFC3 PSF 
sizes are set to 3'.' x 3'.', or the same size as the annulus used in aperture photometry, and the ISAAC and 
NIC3 PSFs are set to 7'.'5 x 7'.'5 to capture any light scattered to high radii. 



3.6. Computing NUV — R Color Gradients for ERS sample 



The selection of the ERS sample of 34 X-ray sources is described in Section 12.41 We select PSFs as 
described in Section 13.51 and use the CLEAN algorithm to deconvolve the images as in Section 13.21 To 
measure color gradients for individual sources, we first compute surface brightnesses in all observed bands 
with the CLEANed images for five concentric elliptical apertures of mean radii 0'.'16, 0'.'33, 0'.'59, 0'.'88, and 
1'.'27. The smallest of these apertures is filled and the larger four are annular. The axis ratio and position 
angle of the ellipses is chosen as described in the first part of Section [31 For each aperture, we interpolate 
through observed B, V, i, z, F125W, and F160W^ surface brightnesses to compute a rest-frame NUV — R 
color. For each source, we convert the apparent angles of arcseconds into physical units of kpc and interpolate 
the radius onto five points: 1.3, 2.67, 4.78, 7.13, and 10.3 kpc. 



4. CHARACTERIZING MEASUREMENT ERROR WITH SIMULATIONS 

Multiple sources of error come into play in the measurements of the colors of AGN host galaxies at 
z ~ 1. First, host colors may be contaminated by the non-stellar colors of AGN emission. In the case of 
integrated colors, photometric measurements with filled apertures directly sum AGN-related emission and 
stellar light. In the case of "outer" colors measured with annular apertures, emission from the central AGN 
may still fall in the photometric aperture due to PSF scattering. Second, our attempted method to correct 
for PSF contamination (described in 13.11) will be sensitive to PSF variation across the field. Third, fiux 
measurements in large annular apertures may be sensitive to sky subtraction errors, depending on the exact 
dimensions of the aperture. 
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These potential errors motivate a detailed approach to estimating measurement errors. Using a suite of 
cloning simulations, we characterize measurement error as function of the strength of a central point source 
and the intrinsic surface brightness of the galaxy host. In these simulations, we modify images of local 
galaxies to approximate their appearances as if they were located at higher redshift. One powerful advantage 
of cloning simulations is that we are able to measure the intrinsic surface brightness of an individual cloned 
galaxy with high S/N before PSF convolution. As this information is difhcult to obtain for high-redshift 
galaxies, cloning will allow us to measure errors at fainter host magnitudes than the alternative approach of 
adding point sources to neighboring galaxies at comparable redshifts to the AGN sample. 



4.1. Implementing Cloning 



Sloan Digital Sky Survey (SDSS, I York et al.l 120001) images of ten local galaxies were downloaded from 



the SDSS website in u'g'r'i'z' filters, mosaiced, and sky subtracted. These ten galaxies were chosen because 
their distribution of inclination angle, half-light radius, and concentration are representative of the "all" 
sample. The natural colors of these galaxies are not representative of the "all" sample, however, so the 
colors are manipulated to match the mean colors of the "all" sample during the cloning process. The final 
reduced image sizes range from 27 kpc to 67 kpc, covering physical regions larger than the diameter of the 
most expansive photometric aperture at high redshift for this sample (25.4 kpc aX z ~ 1.5). The seeing in 
physical units (^ 20 — 50 pc) is several orders of magnitude smaller than the effective image quality expected 
at high redshift (^ 500 — 4000 pc). Ultraviolet images were downloaded from the GALEX public website for 
these four galaxies in both near-UV (NUV) and far-UV (FUV) filters. 

Images of cloned galaxies are generated by inserting scaled copies of these SDSS and GALEX images 
into empty regions of ACS B, V, i, z, ISAAC J, H, Ks, and NIC3 FIWW and F160W images. The change 
in scales from the local universe to high redshift is calculated assuming a flat cosmology with Hq = 70 
km/s/Mpc, i^M = 0.3, and 11a — 0.7 with the IDL "lumdist" function. For an individual cloned instance 
of a galaxy, the target redshift is generated randomly over the range 0.5 < z < 1.5 and used to assign 
GALEX/SDSS filters to the ACS/ISAAC fihers into which they redshift. For a given redshift, the filter 
closest in wavelength is chosen. At the lowest redshift of z = 0.5, the SDSS z' image is far bluer than the 
wavelength required to fill the ISAAC Ks band image (~ 1.5 /im). However, the morphological K-correction 
between rest-frame ^ 0.9 ^m and ^ 1.5 /im is expected to be small. 

Once source filters are assigned to target filters, the scaled copies are convolved with point spread 
functions drawn from the target images. Each of 10 empty target locations in the ACS/ISAAC images is 
chosen nearby a reference star, which is used for convolution. The reference star is not used for simulated 
observation in later steps. The images and PSFs are supersampled by a factor of 3 for this convolution and 
brought back to the target pixel resolution. The ISAAC point spread functions possess considerable noise, 
so 2-D Moffat fits of these PSFs are used for convolutions. The radial profiles of the Moffat fits match the 
radial profiles of the fitted ISAAC PSFs to within 10% at all radii less than 3'.'0 for all 10 target locations. 

For each cloned galaxy, the observed z-band surface brightness in the 0'.'75 < r < 1'.'5 elliptical aperture 
is set randomly between 21 and 31 magnitudes per square arcsecond. The amplitudes of the images in each 
of the other filters are chosen according to the mean rest-frame SED of the actual host galaxy population as 
measured from this data set {{u — i?)" ~ 1.6; {R — J)^ ^ 1.1). The intrinsic host surface brightnesses (i.e., 
before convolution) in the elliptical aperture are computed and recorded for comparison later. 

To test the robustness of wide aperture photometry against contamination from AGN emission, a point 
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source of random amplitude is added to the host galaxy center. The Point Source Fraction, or ratio between 
the flux in the point source and the flux in the entire host galaxy, ranges from 0.1 to 1000. 



4.2. Calculation of Critical Errors 

The images of ten local galaxies are cloned 1600 times in each of 10 insertion locations in the ACS and 
ISAAC fields, 800 with no central point source, and 800 with central simulated point sources. The redshifts 
and bolometric galaxy luminosity vary as described above. For each clone, full aperture photometry is 
performed as described in Section 13.11 using an aperture correction to compute uncontaminated surface 
brightnesses. We compare the measured surface brightness with the true, instrinsic surface brightness in the 
correct aperture before convolution with the PSF. 



4-2.1. Galaxies without Central Point Sources 

Here, we estimate photometric error in the "outer" aperture as a function of intrinsic host surface 
brightness using cloned galaxies that lack central point sources. Given the measured surface brightnesses for 
the "all" sample, we estimate photometric errors and evaluate the reliability of the "outer" photometry for 
this sample. 

Figure [3] plots the difference between the true surface brightness and measured surface brightness as a 
function of true surface brightness for 800 cloned galaxies in all observed bands. The differences are plotted 
in relative flux units, or {fmeas — ftrue)/ fmeas- Shaded regions denote the la random error inferred by fitting 
Gaussian functions to the relative flux difference distributions. The red vertical lines denote the surface 
brightness limit where the implied random photometric error exceeds 0.4 magnitudes. We choose this limit 
for two reasons: First, the expected number of sources with negative flux measurement (a catastrophic 
failure) is limited to a few percent of the subsample whose random errors are of order ^-^40% in flux units. 
Second, 0.4 magnitudes of random error will average down to our expected level of systematic error (^ 0.1 
mag) in samples with ~ 16 — 32 sources. 

Histograms in all panels of Figure [3] show the distribution of measured surface brightnesses in each 
band for the "all" sample. These histograms indicate that only a few sources are fainter than the limit 
corresponding to 0.4 magnitudes of error (red line) in all observed bands redder than B. In the i?-band 
panel, a blue histogram denotes the distribution of B surface brightnesses for sources in the "all" subsample 
with z < 0.638. Sources below this redshift require the i3-band to interpolate a rest-frame u flux; sources 
above this redshift do not require _B-band measurements. Of all the observed bands, the highest photometric 
errors are seen in the _B-band for the "all" sample. However, none of the sources for which i3-band is 
necessary to compute a rest-frame u flux have a photometric error greater than 0.4 magnitudes. 

We conclude that, for the "all" sample as selected, rest-frame u — R color can be reliably constrained 
for all sources that lack central point sources. 



17 



Fig. 3. — Plots of relative flux differences [fmeas — ftme)/ fmeas in ACS and ISAAC bands for cloned galaxies 
without central point sources, plotted as a function of true galaxy surface brightness in the "outer" aperture 
in units of AB magnitudes per square arcsecond. Each filled circle denotes a separate cloned instance of a 
local galaxy. Different colors identify the original local galaxy. Gray dashed lines indicate the mean inferred 
systematic error as a function of surface brightness for these simulations. Shaded regions denote the la 
random error inferred from these simulations. Red vertical lines plot the surface brightness at which the 
total inferred photometric error exceeds 0.4 magnitudes. Solid black lines plot histograms of measured surface 
brightness for the "all" sample in each band. The blue line in the upper left plot shows the histogram of 
measured _B-band surface brightness for galaxies with redshifts z < 0.638, below which _B-band measurements 
are required to compute a rest-frame u flux. Note that very few galaxies exceed 0.4 magnitudes of error in 
any band, except for in the B-band. However, for those galaxies that require an observed B flux to compute 
a rest-frame u flux, none exceed 0.4 magnitudes of photometric error. 
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4-2.2. Galaxies with Central Point Sources 

It is expected that the presence of central point sources will cause error in surface brightness measure- 
ments, although the use of a wide annulus (inner radius 0'.'75) will reduce this effect. As nonthermal AGN 
emission tends to be bluer (rest-frame u — R^ 0.3) compared to galaxy light (rest- frame u~ R^ 1 ~ 2), this 
error is likely to be larger in the observed B and y-band for a sample with 0.5 < z < 1.5. We now quan- 
tify the photometric error in each band as a function of central point source strength through our cloning 
analysis. 

For this analysis, the cloning simulation of section 14.2.11 is repeated, but with artificial point sources 
added of varying amplitudes. The ratio of the total amplitude of the point source to the total ampltitude 
of the galaxy (measured over a 3'.' diameter circle), or Point Source Fraction, is varied between 0.1 and 
1000. Aperture photometry is performed on each of the cloned galaxies and corrected for contamination as 
described in Section [01 We then compare the corrected surface brightnesses to the intrinsic host surface 
brightnesses to evaluate photometric error as a function of Point Source Fraction. 

The relative flux differences in all bands are shown in Figure SI as in Figure [3J but now plotted against 
the logarithm of the Point Source Fraction. Shaded regions denote the la random error inferred by fitting 
Gaussian functions to the relative flux difference distributions. The red vertical lines denote the Point Source 
Fraction at which the random photometric error exceeds 0.4 magnitudes. There are clear trends between the 
reconstructed surface brightness and the amplitude of the central point source in all bands. The photometric 
error is less than 0.4 magnitudes for sources with Point Source Fraction less than ^25 — 50 in each band. 
Notice that the HST V and i bands have the least photometric error for a given Point Source Fraction and 
the B and K bands have the most. 

Figure |4] indicates that "outer" aperture flux measurements will be reliable for sources with Point Source 
Fraction less than 30 in all bands. While the Point Source Fraction is difficult to measure for non-simulated 
data, a more easily- measured proxy can be used. Here we use the ratio between the amplitude in a O'.'IS 
aperture and a circular, annular aperture of inner radius O'.'IS and outer radius 1'.', or Cg.is/i.o- We define 
the Normalized Concentration Co.i8/i.o,norm ^^ 

„ _ ^0.18/1.0, image 

'^0.18/1.0. norm — ~p^ ■ 

L'0.18/1.0,PSF 

This parameter will correct for PSF variation in a field. Figure [5] plots the relationship between 
C'o. 18/1.0, norm and the known Point Source Fraction as measured for the cloned galaxies above in observed 
V . The tightness of the correlation suggests that measurements of Co.is/i.o.norm in actual galaxies can be 
used to look up the Point Source Fraction for a given band. 



4.3. Correction for Contamination in Integrated Fluxes 

Gonstructing GMDs that are corrected for point source contamination requires an unbiased estimate 
of the host galaxy's absolute magnitude. In practice, this is more difficult than estimating uncontaminated 
color in the "outer" aperture, as it requires measurement of the integrated galaxy flux beneath a central 
point source. Here, we describe our method of using imaging to estimate point source contamination for a 
given band. 
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Fig. 4. — Plots of relative flux differences {fmeas — ftrue)/fmeas in ACS and ISAAC bands for cloned galaxies 
with central point sources, plotted as a function of the Point Source Fraction. Each panel plots simulated 
points for a different observed band. All symbols are defined as in Figure [31 
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Fig. 5. — Plot of measured log(Co.i8/i.o, norm) as a function of input Point Source Fraction for the 800 
cloned galaxies in this simulation. Dashed lines mark fitted spline functions. The tightness of the correlation 
suggests that the easily measured Cq.is/i.o is a good proxy for point source contamination. 
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The tightness of the correlation in Figure [5] suggests that Co.is/i.o.norm may be used as a proxy for point 
source contamination. To estimate the point source contamination, we first measure Cp.is/i.o in ^U bands. 
Next, for each band, we fit separate sphne functions to the Co.i8/i.o,norm / Point Source Fraction relations 
in the simulated galaxies (e.g., Figure [S|). These spline functions are inverted to create monotonic look-up 
tables. Finally, the Co.i8/i.o,norm values are used to look up the Point Source Fraction in each band. The 
Point Source Fraction is used to calculate and subtract a correcting factor from the measured integrated 
photometry. 

We estimate rest-frame _R-band absolute magnitudes by performing this procedure in the observed JHK 
bands. These are displayed in the next section in corrected CMDs for the "all" sample. 



4.4. Implications for Observations 

Our cloning studies have constrained systematic and random measurement errors as a function of intrin- 
sic surface brightness and Point Source Fraction in all bands. We use these relations to calculate systematic 
and random errors for each AGN host in the "all" sample. We estimate Point Source Fraction as described 
in Section 14.31 We note that using measured instead of intrinsic surface brightnesses to look up errors will 
underestimate error bars for some of the faintest sources. However, the original sample cut of observed 
-R < 24 (total galaxy magnitude) has largely removed intrinsically faint sources from the sample, and very 
few of the remaining sources have photometric errors larger than 0.4 magnitudes. 



5. RESULTS 

5.1. Corrected Color- Magnitude Diagrams 

Figure [5] plot "outer" u ~ R colors against corrected Mn for the "all" AGN sample. Rest-frame Mfi 
is estimated from observed JHK integrated photometry (S'.'O filled aperture) via interpolation. Corrections 
for point source contamination are measured using source concentration as described in Section 14.31 above. 
Figure [B] shows uncorrected u — R and M/j as light blue and light red points. These uncorrected loci are 
linked with the corrected photometry with vectors. 

The simulations above have shown that photometric measurements using annular apertures avoid color 
contamination due to central point sources reliably. However, estimating uncontaminated absolute magnitude 
is more difficult, as it involves estimating the flux directly beneath the AGN point source. For large point 
source contamination (Point Source Fraction > 30), this method produces uncertain results. Therefore, we 
make no conclusions about the AGN samples involving the mean or individual values of Mb,. 

It is clear in Figure[6]that the difference in mean rest-frame u—R color between soft and hard sources seen 
in Figure [2] is drastically reduced by using annular photometric apertures and correcting for contamination 
from AGN emission. The mean colors are also listed in Table [TJ note that the use of "outer" apertures 
reduces the soft-hard color difference by a factor of nearly three. Given the measurement errors, the mean 
rest-frame u — R color of the "all-soft" and "all- hard" subsamples are statistically equivalent. 

Note that the mean color and magnitude of the obscured population is different in Figure [2] and El despite 
the lack of point sources in this population. The corrected "outer" magnitudes of the obscured population 
are ~ 0.1 bluer than the uncorrected, integrated magnitudes. Because of the negative color gradients in the 



22 



Fig. 6. — Rest frame u — R vs. Mji color-magnitude diagram, corrected for point source contamination. 
Symbols are plotted as in Figure [H u — R colors have been measured with annular "outer" photometric 
apertures. Mji magnitudes are measured with integrated apertures. Sources judged to be marginally con- 
taminated by a central point source associated with AGN emission (Point Source Fraction in observed B 
greater than 0.35) are corrected for this contamination in Af^. For these contaminated sources, uncorrected, 
integrated measurements are shown as faint, empty points. Uncorrected, integrated loci are connected to 
corrected, "outer" loci with dashed lines. Rest-frame values are estimated from observed photometry as in 
the text. Blue and red error bars denote the mean colors and luminosities for the soft and hard sources, 
respectively, with errors calculated as in Section 15.41 Small black circles denote the colors and magnitudes 
of the underlying galaxy population, selected as explained in Section 12.4.11 Note that photometry is mea- 
sured with 3" filled apertures for this control sample, not "outer" apertures. The soft source at the bluest 
extreme of the diagram (XID 100) has a significantly bluer "outer" color than its integrated color, unlike the 
remainder of the "all" sample; it appears to be composed of a blue point source with a bluer star-forming 
extension to the NE. 
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host galaxies, the outer regions are bluer than the nuclear regions. For the unobscured sources, this effect is 
overwhelmed by the removal of the central point source contamination, which causes the corrected "outer" 
magnitudes to be redder than the integrated magnitudes. 



5.2. Mean SED Comparison 

The surface brightnesses of the AGN sample are plotted as de-redshifted spectral energy distributions 
(SED) in Figure [T] Surface brightnesses are plotted in absolute AB magnitudes, modified to account for a 
de-projection of inclined disks (i.e., fluxes are multiplied by 5/a, where b is the semiminor axis and a is the 
semimajor axis) . Figure [7] plots the logarithmic mean (or geometric mean o f magnitudes) of these SEDs for 
the "all-soft" and "all-hard" subsamples. SEDs of lBruzual fc Chariot! |2003l) Single Stellar Population (SSP) 



bursts are shown for comparison. The SEDs of the soft and hard sources appear to be similar, although the 
soft sources are ^ 0.3 magnitudes brighter than the hard sources at all wavelengths. 

The mean SEDs of the two populations are inconsistent with a starburst population of age ~ 25 Myr. 
The SEDs more resemble those of intermediate and old stellar populations, although the rise at near-infrared 
wavelengths compared to these SSP models implies that dust extinction is present. Young stellar populations 
and their accompanying ultraviolet emission may be hidden by this dust. 



5.3. Mean Color Gradients in ERS sample 

The principal advantage of using the "outer" aperture is the avoidance of central point sources; the 
primary disadvantage is that is does not probe the central regions of host galaxies. The inner radius of the 
"outer" aperture was set by the image quality of the ISAAC JHK imaging, approximately 0'.'65 in the worst 
case. We now utihze F125W and F160W imaging in GOODS-South from the WFC3/IR imager aboard 
HST, with a spatial resolution of ^ 0'.'2, to constrain host colors at smaller radii. Although the spatial 
coverage of the WFC3 infrared imaging is less than 50 square arcminutes, the improved depth relative to 
ISAAC {Hab — 27.0 for point sources, 5cr) permits extending the sample to fainter magnitudes. 

We compute NUV — R color gradients for each host galaxy as described in Section [3T6l Figure [8] displays 
the mean NUV — R color gradients for the "ERS-soft" and "ERS-hard" subsamples. The shaded regions 
denote the 68% confidence intervals as computed in Section \5A\ below. Two curves are displayed for soft 
sources; the blue curve shows the mean color gradient for all 15 soft sources and the hatched curve shows 
the mean color gradient for the 10 soft sources whose Point Source Fraction in V-band is less than 10. Point 
Source Fraction is measured in the y-band as described in Section [4.31 



Note that because the images are CLEANed and the central aperture is filled, contaminating point 
sources will affect the color of the innermost point only (r = 1.3 kpc) in Figure [H The innermost point 
of the mean color gradient of the full soft sample, denoted by the blue curve, is significantly bluer than all 
other points. However, when sources with significant point source contamination (Point Source Fraction > 
10) are removed from the soft sample, as shown with the hatched curve, the innermost point returns to a 
color that is similar to the colors at all other radii. The removal of these contaminated sources does not 
affect the mean colors of the hosts beyond r ~ 2 kpc. 

The colors of the hard and soft sources are statistically identical at all radii beyond r ~ 2 kpc whether 
objects with significant point source contamination have been removed or not. Removing objects with 
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Fig. 7. — Logarithmic mean of surface brightnesses for all-soft sample and all-hard sample. The dashed lines 
represent the one-sigma uncertainties on the mean SEDs. Red lines denote measurements for soft AGN and 
blue lines denote measurements for hard AGN. Error bars are measured from full cloning simulatio ns and 
a bootstrap procedure to estimate error due to finite sampling. SEDs of iBruzual &: CharlotI ([20031) Single 
Stellar Population (SSP) bursts of varying ages are shown for comparison. 
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Fig. 8. — Mean color gradients for three populations in the WFC3 ERS sample: Obscured sources (HR 
> —0.2), all unobscured sources (HR < —0.2), and only unobscured sources with little point source con- 
tamination (V Point Source Fraction < 10). Filled regions are 68% confidence intervals. With the removal 
of sources strongly contaminated by central point sources, the mean color gradients of the obscured and 
unobscured sources are similar to ^ 0.5 magnitudes at all radii greater than 1 kpc. 
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central point sources from the soft sample allows us to make an unbiased color comparison between soft and 
hard sources for the central r = 1.3 kpc point; when these objects are removed, the remaining soft sample 
and the hard sample are statistically equivalent in color at r = 1.3 kpc. Assuming that the point source 
contamination in the objects with Point Source Fraction > 10 is caused by nonstellar emission associated 
with the AGN, we may conclude that the colors of the soft and hard sources for this subsample are similar 
at all radii greater than r ^ 1 kpc. 



5.4. Error Considerations 

The method of averaging can possibly affect the mean photometric properties of a sample. When the 
averaging is performed in flux space, with equal weights on each photometric point, the principal outcomes 
of the comparisons between the various types of AGN in this sample are unchanged. Similarly, when 
the averaging is performed in flux space with stronger weights towards fainter fluxes, the comparison is 
unchanged. 

Systematic measurement errors are seen in the cloning simulation; we estimate these for individual 
galaxies and subtract them. The random errors on the mean SEDs at a particular wavelength and on the 
mean colors of subsamples are the sum of two components: The first is a random error caused by sky- 
subtraction errors, photon error in the photometric measurement, and other errors associated with basic 
aperture photometry. The second major error, dominant over photometric error, is due to small-number 
st atistics in these su b-samples. We estimate this error with a Monte Carlo bootstrap method, as outlined 
in lPress et al.l (|l992l ). In this method, new sub-samples are repeatedly drawn (with replacement) from the 
overall data set, analyzed, and used to populate a distribution. The random error due to the finite size of the 
sample is then given by the standard deviation of critical parameters (e.g., mean rest-frame u magnitude) 
in this distribution. 



The result of this analysis is that the bootstrap error in u — R is 0.11 magnitudes for the all-soft 
subsample and 0.10 magnitudes for the all-hard subsample. This can be understood in terms of the intrinsic 
brightness distribution and color scatter in the samples. The intrinsic variation inu — R color for this sample 
is (T '-^ 0.7 magnitudes. The error on the mean of some parameter due to incomplete sampling of a population 
goes as a/\/N, where a is the standard deviation of the parameter for the parent population and N is the 
size of the sample. With A'' ~ 31 — 42 for the all-soft and all-hard subsamples, respectively, the error on the 
mean is '^ 0.13 and ^0.11 magnitudes for u — R color, respectively. These alternate estimates are nearly 
the same as that estimated with a Monte Carlo bootstrap procedure. 



5.5. Are Obscured AGN Hosts More Dusty than Unobscured AGN Hosts? 



u — R and NUV — R colors are sensitive probes of low levels of star formation (jSalim et al.ll2007[ ). but 
these colors can be affected by dust extinction. Here we search for systematic differences in galactic dust 
extinction between the obscured and unobscured AGN host populations. In a rest-frame U — V , V — J 
color-color diagram, the effects of dust extinction can be separated to some extent fr om the effects of st ellar 



age (see Figure |9]). V — J color is more sensitive to dust extinction than stellar age (jWuvts et al.ll2007l ). 



We compute V — J "outer" colors to estimate the differential effects of dust extinction. Rest-frame J 
band magnitudes are calculated from observed BVizJHK photometry using linear interpolations as described 
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in Section[3l Half of the sample has z > 0.85, beyond which rest-frame J shifts redder than observed Ks, the 
reddest band available in this sample. Color extrapolations are required beyond this redshift. We estimate 
the error in this extrapolation by measuring the natural dispersion in rest- frame z— J "outer" color (observed 
H — Ks) at a single rest-frame R — z (observed J — H) color slice for sources with z < 0.85. This dispersion 
is 0.10 magnitudes, implying that the extrapolation error is less than 0.1 magnitudes for all sources in the 
sample. The error on the mean of the samples will be much smaller than this. 

The bottom panel of Figure IH] plots the "outer" colors of the "all" sample in a rest- frame U —V, V — J 
color-color diagram. Only sources for which observed Ks is measured are included, i.e., only sources with 
ISAAC imaging in the southern GOODS field. The mean rest-frame U —V color of the two populations are 
discrepant by 0.16 magnitudes, which is consistent with the results shown in Table [1] within the errors. The 
mean V — J color of the obscured population is redder than the unobscured population by 0.17 magnitudes. 
These differences are not statistically significant, suggesting that the two populations are similar in mean 
galactic dust extinction. The two populations do not occupy clearly separated regions in the diagram. 
Approximately one- half of the sample falls in a region corresponding to Ay > 1, suggesting that a significant 
fraction of the sources a re dusty, which is consistent with other studies using similar data and methodology 
(jCardamone et al.ll2010( ). 



5.6. u — R Properties of AGN Hosts 

As seen in Figure [7] and Table [1] the mean SEDs and mean rest-frame u — R colors of AGN hosts beyond 
r ^ & kpc are largely similar across X-ray obscuration. The hosts of the soft AGN, on average, are ~ 0.11 
magnitudes bluer \nu — R and ~ 0.2 magnitudes bluer in NUV — R than the hosts of the hard AGN beyond 
r '-^ 6 kpc. This effect is weak and not statistically significant. When this analysis is repeated for luminous 
objects with Lx > 10^"^ ergs s~^, the u — R color of the soft sources is statistically equivalent to that of 
the hard sources. The soft sources are 0.17 magnitudes bluer than the hard sources, but this effect is not 
statistically significant. 

We have checked this result with a simple stacking exercise, in which aperture photometry is performed 
on a stack of the processed images of a particular population. For this procedure, we use sky-subtracted 
images whose ACS/ISAAC alignment is good to ^ O'.'l that have been deconvolved with CLEAN as described 
in section 13.21 This stacking method simulates the measurement of fiux-averaged colors, which effectively 
weights brighter sources higher than fainter sources relative to a logarithmic mean. We use observed V and 
J colors, which correspond to nearly u and R aX z ^ 0.8. Using a single circular, annular aperture of 0'.'75 
inner radius and 1'.'5 outer radius, the observed Vab — Jab colors of the stacked soft and hard samples are 
0.81 and 0.95 magnitudes, respectively, suggesting that the soft sources are 0.14 magnitudes bluer than the 
hard sources beyond r ~ 6 kpc. Although this value is larger than the difference in u — R colors measured 
via SED averaging, the similarity is reasonable considering that (1) the bands are in the observed frame and 
not the rest-frame and (2) a single aperture is being used for a variety of host morphologies. 

We have seen that the mean rest-frame V — J color is also similar for the soft and hard sources (Figure 
[9|). This implies that the mean rest-frame u — R colors of the soft and hard subsamples are affected by 
dust extinction at similar levels. Given both the similarity of the rest-frame u — R colors and the similarity 
of the mean dust extinction in the soft and hard subsamples, it follows that the mean extinction-corrected 
rest-frame u — R colors of the soft and hard sources would also be similar. 
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Fig. 9. — Top panels: Rest-frame U — V vs. V — J color-colo r diagrams of all ga laxies in a mid-IR selected 
sample with Ly > 5 x lO^i© (included with permission from lWuvts et al.l (|2007l )). SDSS+2MASS galaxies 
(small gray dots) are plotted as a local reference. Top left panel: Galaxies are color-coded by dust extinction. 
Dust extinction affects rest-frame V — J color more than U — V. Top right panel: Galaxies are color-coded 
by mean stellar age. Stellar age affects rest-frame U — V color more than V ~- J. Bottom panel: Rest-frame 
U — V vs. V — J color-color diagram using "outer" aperture for AGN with observed Ks available. Small 
black circles denote the colors of the underlying galaxy population, selected as explained in Section 12.4.11 
Note that photometry is measured with 3'.' filled apertures for this control sample, not "outer" apertures. 
The mean colors of the AGN sample are shown in the center of the diagram; the obscured locus is close to 
the unobscured locus in T^ — J, suggesting that the populations have similar mean extinctions. 



2.0 
1.5 

^ 1.0 

I 

0.5 

0.0 



■Lfij- 



+ 



. :»'■■* ■• 






A„ =0-0 0.5 1.0 1-5 2.0 2.5 



1 2 






■ 1 ■ ■ ■ ■ 1 ■ ■ ■ 


■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ 


2.0 


■ • log(l^ > 43 
r • log(l^ < 43 

■ il HR > -0.2 


-y-: ■'■■■.■^ifc ■ ■ 

..-■■^i:#:!v.*^ 

. 1 .... 1 .... 1 .... 1 .... 1 . ' 


1.5 

1 

r 1.0 
I 

0.5 
0.0 


- 


• 



-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



29 



6. DISCUSSION AND CONCLUSION 



6.1. Comparison with Previous Work 

This work is the first to measure the UV/optical colors of a large sample of AGN host galaxies at 2; ~ 1 
while quantitatively, reliably correcting for point source contamination in individual galaxies. Ignoring or 
insufficiently correcting for point source contamination results in systematic error in comparing the colors of 
hard and soft X-ray sources (see Figure [5]) . 

Direct comparisons of our results with those obtained for Seyfert galaxies in the local universe are 
problematic, as numerous AGN selection techniques are commonly used in the literature. At low redshift, 
the effects of point source contar aination on host galaxy color measurements are negligible because the 
galaxies are well- resolved. Locally. ISchawinski et al.l ( 20091 ) find little host color difference between obscured 



and unobscured X-ray-selected AGN, in agreement with our results. Although the sample is small, their 
selection techniques and use of X-ray-derived obscuration measurements are more analogous to those utilized 
in this study. Given the similarity of our selection criteria, the agreement of our results suggests that the 
mean color offset between soft and hard sources has not evolved since z ^ 1. 



Ho et al.l (|2008l ) find that unobscured broad-line AGN hosts at low-redshift are observed to possess 
large reservoirs of HI gas. This result is interpreted as being evidence agai nst a feedback scenario, in which 
unobscured AGNs are thought to be caught in the act of expelling gas (JHo et al.ll2008l ). No comparison 
is made between Type I and II sources. Considering the use of different selection methodologies and mea- 
surements, there is no i nconsistency betw een this study and our own. For a sample of mid-infrared Type II 
QSOs at 0.3 < z < O.S. lLacv et al.l (|2007l ) find that the host galaxies exhibit significant star formation rates 
(3 — 90 Mq yr~^) and that disk inclination correlates with silicate features, implying that at least some of 
the reddening arises from the host galaxy. However, it is not clear how these star formation rates would 
compare to similarly selected Type I QSOs at similar redshift. Moreover, the AGN selection techniques, 
selected luminosity range, and redshift range are different from those in the present study, preventing direct 
comparison. 



Page et al.l (|2004 ) compares the 850 fim fluxes of luminous {Lx ^ W^^ ergs s~^). X-ray-selected AGN 
at 1 < z < 3, finding that obscured sources have significantly more submillimeter emission than unobscured 
sources. As a function of redshift, the mean flux difference between obscured and unobscured sources is 
consistent with ze ro at redshifts below z ~ 1.5 and rises significantly at higher redshifts. Contrastingly, 
Shao et al.l ( 20101 ) measures the far-infrared-derived star formation rates of X-ray selected intermediate- 



luminosity AGNs {Lx ~ 10'*^ ergs s~^) at 2; ^ 1 in GOODS-North, finding no dependence of star formation 
rate on X-ray absorbing column density. It appears that the mean star formation rate of obscured sources 
diverges from that of unobscured sources at either the highest X-ray luminosities {Lx ^ 10*^^ ergs s^^), at 
higher redshifts beyond z ~ 2, or in both of these regimes. 

CPIO measures the rest-frarne U-V colors of intermediate-luminosity. X-ray-selected AGN in the Ex- 
tended Groth Strip ( Pierce et al.ll2010b[) . They find that soft sources are systematically bluer than hard 
sources in both nuclear (r < 0'.'2), extended (0'.'2 < r < I'.'O), and integrated apertures. The mean difference 
in extended color between soft and hard populations in our GOODS sample differs from the CPIO value 
by 2.5cr (0.3 mag in u — R). Our color offsets are consistent to within la (0.1 mag in u — i?) when the 
CPIO sample is trimmed of all sources that they identify visually as having "definite" or "possible" point 
sources. In addition, color offsets are similar when no correction for point source contamination is used in 
our sample (0.3 mag in u — R; see Figure [2] and Table |T]). In CPlO's Figure lid it is clear that sources 



30 



visually identified as possessing "definite" or "possible" point sources display anomalous color gradients (i.e., 
blue nuclear regions). These lines of evidence point to some amount of point source contamination in CPlO's 
extended colors. 



Cardamone et al.l (|201Cll ) find that X-ray selected AGN hosts at 



1 are frequently in dust-enshrouded 
galaxies. Approximately one-half of our sample occupy a region in the rest-frame U — V vs. V — J diagram 
corresponding to Ay > 1, so our results are not inconsistent. Although the authors do not explicitly compare 
the locations of obscured and unobscured sources in this plot, it is clear from their Figure 2 that the mean 
locations of these populations are similar, as in our sample. 



6.2. Implications for AGN Unification and Galaxy Evolution Scenarios 

In this paper, we have assembled an X-ray-selected AGN sample and measured the mean u — R colors 
beyond r ^ 6 kpc in the host galaxies. We found that the mean rest-frame SEDs were similar for soft and 
hard sources both for objects of intermediate X-ray luminosities {Lx > W^^ ergs s""'^) and objects of higher 
luminosities {Lx > lO'*'^'^ ergs s^^, (Lx) ^ 10^^'^ ergs s^^). We also found that the mean loci of the soft 
and hard populations in a rest-frame U — V vs. V — J were similar, implying that the mean properties of 
galactic dust extinction of the two populations were similar. In Section [5.31 we used high spatial resolution 
WFG3/IR imaging to compute the mean rest-frame NUV ~ R color gradients of soft and hard sources. The 
CLEAN deconvolution algorithm was used to remove PSF contamination effects before computing colors. It 
was seen that the colors of soft and hard sources were significantly different only in the nucleus (r < 1.3 kpc), 
and when sources with strong point sources in the observed y-band were removed from the soft sample, the 
mean colors of the soft and hard sources were similar with respect to the error bars at all radii. The limiting 
spatial resolution of the WFC3 imaging prevents us from probing colors at radii smaller than r '^ 1 kpc. 
We conclude that the host colors of soft and hard sources are statistically equivalent at all radii larger than 
r ~ 1 kpc. These observations indicate that unobscured AGN are not redder than obscured AGN at z ^ 1 
at intermediate X-ray luminosities. 

In our sample, galaxy-wide dust extinction is weakly or not correlated with nuclear AGN obscuration as 
probed by X-ray hardness ratio. Taken with the evidence that galaxy color is weakly or not correlated with X- 
ray hardness ratio, this implies that the conditions of star formation and dust extinction are uncorrelated with 
the condit ions of neutral hydrogen obscurat ion on nuclear scales. These observations favor AGN unification 
scenarios ( AntonuccilllQQStlUrrv et al.lll995l ). in which AGN obscuration is determined by the orientation of 
a torus on parsec scales with respect to the observer. This conclusion appears to be independent of X-ray 
luminosity in our sample for luminosities lower than lO'*'*'^ ergs s~^. 

Recently, several models of galaxy form ation have been able t o explain the bimodality of galaxy col 



ors through AGN feedback processes (e.g., iDi Matteo et al. 
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Lapi et al.l 120061 : iHopkins et al. 
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Menci et al.l 120081 ). These models imply that unobscured AGN would be associated with a quenching of 
star formation, resulting in a reddening of th e stellar population, even in the outer r egions of the host (be- 
yond r ^ 6 kpc) in the most energetic cases ( Hopkins et al.ll2006t iMenci et al.l 120081 ). The color of a single 
stellar population (SSP) reddens by ~ 0.7 magnitudes in u — i? ov er a quasar lifetim e following the ceasing of 
star formation (solar metallicity from age 10 Myr to age 100 Mvr. iMarastonI (J20051 )). and thus the difference 
in color resulting from quenching in unobscured hosts should have been detectable in this study. 

Our observations instead suggest that the majority of intermediate luminosity AGN at z ^ 1 are 
not undergoing nor have recently experienced rapid, galaxy-wide quenching due to AGN-driven feedback 
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processes. A second possible interpretation is that if AGN-related blowout events have occurred in the 
population of intermediate luminosity AGNs at 2 ~ 1, the radial extent of quenched stellar populations 
must be restricted to ^ 1 kpc on average or the color signature of the quenched stellar population must have 
since been rewritten by fresh gas accretion and resulting star formation. 
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